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The influence of the tail features of the local magnetic field probability density function (PDF) 
on the ferromagnetic Ising model is studied in the limit of infinite range interactions. Specifically, 
we assign a quenched random field whose value is in accordance with a generic distribution that 
bears platykurtic and leptokurtic distributions depending on a single parameter r < 3 to each site. 
For r < 5/3, such distributions, which are basically Students and r-distribution extended for all 
plausible real degrees of freedom, present a finite standard deviation, if not the distribution has got 
the same asymptotic power-law behavior as a a-stable Levy distribution with a = (3 — t)/(t — 1). 
For every value of r, at specific temperature and width of the distribution, the system undergoes 
CO ■ a continuous phase transition. Strikingly, we impart the emergence of an inflexion point in the 

temperature-PDF width phase diagrams for distributions broader than the Cauchy-Lorentz (r = 2) 
which is accompanied with a divergent free energy per spin (at zero temperature). 

o 



I. INTRODUCTION 



Disorder is ubiquitous in Nature. Regarding materials and their statistical properties, disordered magnetic systems 
have been systematically studied in condensed matter and statistical physics. From a theoretical point of view, the 
most studied case has certainly been the Random Field Ising Model (RFIM) fH, because of its simplicity as a 
frustrated system and relevancy to experiments [H, H[ which has been quite boosted after the identification of the 
*^ . RFIM with diluted antiferromagnets in the presence of a uniform magnetic field [B, [B, [B, 0] and several ferromagnetic 
' compounds as well [B, 0, 11] • 

In order to generate the local random field, both the Gaussian and the bimodal probability density function (PDF) 
1 ^ 1 ■ have intensively been used [B, [ToL [TTj] . Nevertheless, controversy over the order of the low-temperature phase transition 
has still been at the helm of several discussions. On one hand, a high temperature series expansion up to 15th order 
showed a continuous phase transition for both the Gaussian and the bimodal PDF TBl- On the other hand, from an 
exact determination of the ground states in higher dimensions (d = 4), Swift et al [H found a discontinuous phase 
transition for the bimodal random field, whereas for d = 3 dimensions and the Gaussian distribution the transition is 
continuous. By applying the Wang-Landau algorithm [ldj], recent simulations on 3D lattices claimed the discovery of 
ly-^ , first-order-like features in the strongly disordered regime for both those PDFs [IB, EB] • 

As an alternative to the above mentioned approaches, there is the mean field theory which can present a good 
qualitative agreement with some short-range interaction models and experiments. Once more, the Gaussian and the 
bimodal PDF have been widely investigated [13, fT§l as well as related distributions such as the trimodal [IB, HB] and 
the double-Gaussian [H| or the treble-Gaussian [23| . In the Gaussian RFIM case, the phase diagram only presents 
continuous phase transitions , whereas in the bimodal case the phase diagram presents a continuous phase transition 
for high temperatures and low random-field intensities and for low temperatures and high random-field intensities a 
first-order transition arises therefrom [HI]. In other more elaborated cases a rich critical behavior can be found for 
finite temperatures as it has been recently conveyed in [HI, [IB] . Accordingly, we can understand that the choice of the 
local random field PDF is of crucial importance for a good theoretical description of real systems. In this particular 
context and based on the identification of the RFIM with diluted antiferromagnets in a uniform field, for which the 
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local random fields are expressed in terms of quantities that vary in both signal and magnitude [5j, the use of 
continuous PDFs has demonstrated to be a very promising approach [2~fl , [HJ . 

The utilization of Dirac Delta and Gaussian related distributions is much supported on the easiness of the analytical 
treatment of the subsequent equations as well as the pervasiveness of the Gaussian distribution. Although the Gaus- 
sian was assumed for many generations as the "natural distribution" , in the last decades the concept of (asymptotic) 
scale- invariance of probability density functions have abundantly emerged [231 ] . In the realm of disordered systems, 
PDFs different to the n-Gaussian or the n-Dirac Delta were used to explain the critical behavior of several com- 
pounds. For instance, PDFs with very fat tails were introduced to analyze organic charge-transfer compounds like: 
N-methyl-phenazium tetra-cyanoquinodimethanide (NMP-TCNQ), quinolinium-(TCNQ)2, acridinium-(TCNQ)2 and 
phenazine-TCNQ, as first reported in Refs. (24[. Conversely, a sub-Gaussian distribution was used to account for the 
magnetic properties and the critical behavior of poly(metal phosphinates) [2a ]. Last but not least, as was proven 
by Gosset [26|, asymptotic scale invariant distributions can be derived from the Gaussian distribution when finite 
elements are taken into account so that finite and scale-dependent systems can be treated as infinite and (asymp- 
totically) scale- independent. Therefore, the study of more general continuous PDFs turns up very interesting as it 
furnishes a more widespread picture of disordered magnetic systems than the distributions used up to now. With such 
a goal in mind, we study herein the aftermath of applying a more general family of continuous PDFs in the mean field 
RFIM. Explicitly, our PDF reproduces the r- and i-distributions for real degrees of freedom. For specific values of 
the triplet composed of the degree of freedom, the temperature and the PDF width, our results show that the system 
experiences a continuous phase transition that does not dependent on the finiteness of the standard deviation and the 
scale behavior (dependence or independence) of the random field. Moreover, for PDFs fatter than the Cauchy-Lorentz, 
we determine the emergence of an inflexion point in the temperature versus PDF width phase diagrams that coexists 
with a divergence at zero temperature of the free energy per spin. 



II. THE MODEL 

The infinite-range-interaction Ising model in the presence of an external random magnetic field is defined in terms 
of the Hamiltonian, 
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where the sum Y2uj) runs over a U distinct pairs of spins Si = ±1 (i = 1,2, N). The random fields {Hi} are 
quenched variables and ruled by a PDF that is defined by a parameter r (generic degree of freedom). For r < 1, 
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(with \H\ < [B T (1 — r)] 1 ^ 2 ) which is the generalized r-distribution, and for r > 1, we have 
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which is the generalized Student-t distribution. By generalized we mean that the degrees of freedom, to and n, of t- and 
r-distributions are extended to the entire domain of feasible real values according to the relations t = (to + 3) / (to + 1) 
[to > 0] and t — (n — 4) / (n — 2) [n > 2], respectively. In Eqs. and ([3]), T[.} is the Gamma function and B T is 
given by 



B T = 
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where u> is the width of the PDF. For r < 5/3 the width and the standard deviation, <j, are related by 

(5-3r)o- 2 = (3-t)lu 2 . (5) 



Alternatively, the functional form of Eqs. ([2|) and ([3]) can be obtained by optimizing the entropic form presented in 
[27l ] by applying the concept of escort distribution, p (H) = P T (H) / J P T (H) dH [28|, [29| , and for that is many times 




FIG. 1: (Color online) Random-field probability distributions for some values of the parameter r (from botton to top: r = 
5/2, 3/2, 1, 1/2 and —20), in the normal (left panel) and log-log scale (right panel). We have used lu — 1 in all cases. 

called q-Gaussian. In this case, oj 2 plays the role of the constraint, J H 2 p(H) dH — ui 2 , which is always finite for 
r < 3 with the corresponding Lagrange multiplier given by Eq. (|4|). Expressly, u> 2 represents the standard deviation 
of the escort distribution and it is finite even when the distribution per se has got a divergent standard deviation, 
J H 2 P(H) dH — a 2 . Therefore, it represents a way of appraizing the broadness of the distribution and this is the 
reason why we named lj width. Recently, Pi( s ) (H) has also been coined generalized Lorentzian [301 ] . Although we 
acknowledge both nomenclatures we use the traditional terminology of r— and t— distributions that is quite well 
established in the Statistics community since a long time. The PDF defined in Eqs. @ and (J3]) is symmetrical 
around H = and represents a family of continuous distributions that recovers some well-known distributions using 
appropriate limits, namely: 

• the uniform distribution, for r — > — oo; 

• compact support distributions (limited), for r < 1; 

• the Gaussian distribution, for r — > 1; 

• the Cauchy-Lorentz distribution, for r = 2. 

• Dirac Delta, for every r < 3 and u> — > 0. 

To boot, the functional form ^ is an asymptotic power-law decaying PDF with finite standard deviation for 
1 < t < 5/3 and an asymptotic power-law decaying PDF, but with infinite standard deviation instead. In both cases 
the decay exponent is equal to 2/(r — 1). The latter case is also capable of reproducing the tail behavior of a-stable 
Levy distributions 

/oo 
exp[-a \k\ a +ikH] dk, 
-oo 

with a = (3 — t) / (r — 1) and broadness a, whose escort-distribution has got a finite width as well. For the case of 
the Cauchy-Lorentz, a = 1 (t = 2) [the only case for which Levy distributions are explicitly defined in real space], 
the parameter a is equal to width u>. Accordingly, if we bear in mind the previous work by Aharony we can 
hold that our enquiry also sheds light on the low temperature behavior of the random-field Ising model with the local 
magnetic field associated with a a-stable Levy distribution. In Fig. [TJ we depict PDFs © and §S§ for some values of 
r. Regarding the kurtosis, 
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the distribution is platykurtic, k < 3, for r < 1 or leptokurtic, k > 3, for r > 1. At this point it is important to stress 
that, as it has been made until now, in spite of being able to present non-mesokurtic distributions the combination 
of Gaussians results in asymptotic scale-dependent distributions. 

From the free energy, F({Hi}), associated with a given realization of site fields, {Hi}, we calculate the quenched 
average, [F({Hi})] H , 

The general mean field result of the free energy per spin, in terms of any PDF of the random fields, is well-known 
[l7L [la ] , and is given by 

/-ym 2 -^( log[2 cosh 8{ Jm + H)]) H (8) 

and the magnetization is given by, 

m = {twh\3(Jm + H)]) H , (9) 
where ( . . .)h stands for averages over realizations of the disorder, i.e., 
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Close to a continuous transition between ordered and disordered phases, the magnetization m is small. So, we can 
expand Eq. ^ in powers of m, 



where the coefficients are given by 



with 



m = Am + Bm 3 + Cm 5 +0(m 7 ), (10) 



A = /3J{l- Pl }, (11) 

B = - { -^{l-A Pl + 3p 2 }, (12) 
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C = ^H2- 17^+30,92-15,93}, (13) 



Pk = (tanh 2fc {pH)) H 



With the aim of finding the continuous critical frontier we set A — 1, provided that B < 0. If a first-order critical 
frontier also occurs, the continuous line must end when B = 0; in such cases, the continuous and the first-order critical 
frontiers converge at a tricritical point, whose coordinates are obtained by solving the equations A = 1 and B = 0, 
on condition that C < 0. Thus, for A = 1, we obtain 

kT 

— = l-(tanh 2 (/?#))„ . (14) 

In the following section, we discuss the role of PDFs ((2]) and (j3]) when they are considered in the formulae presented 
in this section. Our survey includes the analysis of the phase diagrams for the whole domain of r. 



III. FINITE TEMPERATURE ANALYSIS 



Following the above presented results, we proceed by calculating the critical frontiers of the model when the 
temperature is different from zero. In the RFIM, we have a single transition between the two possible phases of the 
magnetization: the ferromagnetic phase (m ^ 0) and the paramagnetic phase (m = 0). The critical frontier separating 
these two phases is found by solving Eq. (fT4"|) . On account of the fact that Eq. (fT4"|) is analytically unsolvable, we have 
been compelled to solve it by numerical means using the Global Adaptative Strategy algorithm 13 II that has been 
proven as the best (i.e., fast and accurate) numerical integration procedure for smooth integrands [321 ]. 
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FIG. 2: (Color online) Phase diagram of the model, in the plane temperature vs to (in units of J), for some values of the 
parameter r < 1. The grey dotted line is for r = — oo and u> (T = 0) = J; the brown dashed line is for r = —20 and 
lo (T = 0) = 0.9831 ... J; the dot-dashed orange line is for t — and lu (T = 0) = 0.8660 ... J; the dotted green line is for 
r = 1/2 and lu (T = 0) = 3-\/5/64J and the black full line is the Gaussian case with oj (T = 0) = y/2/-!rJ. We can observe 
continuous phase transitions between the Ferromagnetic (F) and the Paramagnetic (P) phases for all values of r. The points 
onto the u>/J axis were exactly calculated through a zero-temperature analysis [section IV. A] where from we can see a good 
agreement between the analytical and numerical results which by interpolation indicates discrepancies never greater than 1%. 



A. Platykurtic case: r < 1 

Let us denote /j and m, as the free energy and the magnetization for this regime of r, respectively. Thus, Eqs. ([5]) 
and © become, 

t2 i c — 1 — 
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and 



m, = / %/Bt(1 t> Pi(if) tanh/3(jTO, + ff) , (16) 

where Pi{H) is given by the PDF in Eq. |J2J|. The continuous critical frontier has been found when we have solved 
Eq. (|14[) . For all solutions obtained, we have calculated a negative value of B, Eq. (fT!?|) , which has confirmed the 
continuous character of the phase transition. 

If a first-order transition existed as well, the critical frontier would be found by equalizing the free energy at each 
side of this line, i.e., f(m = 0) = /(to =/= 0). Using this procedure, we have numerically determined the critical 
frontiers separating the paramagnetic and ferromagnetic phases, for typical values of r < 1. We have confirmed that 
the above coefficient B, Eq. (fP2|) . is always negative. The phase diagram is shown in Fig. [21 on the plane defined 
by the temperature, T, and the PDF width, to (both in units of J), for some typical values of r < 1. In that 
figure, the lines represent the numerical solution of Eq. (|14|). whereas the points were analytically obtained through 
a zero-temperature analysis, which is going to be discussed in the next section. Notice that the ferromagnetic phase 
is reduced by increasing the parameter r from r = — oo to r = 1 as shown in Fig. and for the maximum value for 
r-distributions, r = 1, we recover the simple phase diagram of the Gaussian distribution [13]. 
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FIG. 3: (Color online) Phase diagram of the model, in the plane temperature versus uj (in units of J), for some values of 
the parameter r > 1. The black full line is the Gaussian case with lo (T = 0) = \f^J\ red dashed line is for r = 3/2 and 
w (T = 0) = "^f - ^' t ne purple dotted line is for t — 2 and u (T = 0) = ^ J; the dot-dashed blue line is for r = 5/2 with 
oj (T = 0) = 0.4754 ... J. We can observe continuous phase transitions between the Ferromagnetic (F) and the Paramagnetic 
(P) phases for all values of r. The points represent the results obtained by the zero-temperature analysis. Notice the change 
in the concavity of the critical frontier for large values of r (> 2.0). The vertical dashed line is u) = 0.275J which is close to 
the inflexion point of the critical line for r = 5/2. In this figure, we have distinguished the points with finite free energy per 
spin from the points with a divergent free energy per spin representing the latter by empty circles. Again, we can observe a 
good agreement between the numerics and the expansion at T = 0. The difference between the analytical approximation and 
interpolation is again never greater than 1%. 



Analogously to the platykurtic case, we denote f s and m s as the free energy and the magnetization per spin for 
this regime of r. The expansion of the magnetization Eq. (fTU|) is valid for this case as well, but the averages over the 
disorder (...)# must be made according to PDF fl3j, 



where in this case the integration limits are taken in the range (-co, +oo). 

By considering PDF ([3]), the above presented procedure for the determination of the critical frontiers can be 
employed once more. In other words, Eq. (|14p provide the continuous critical line of the phase diagram. Using this 
procedure, we have numerically evaluated the critical frontiers separating the paramagnetic and ferromagnetic phases 
for typical values of r > 1. Like the platykurtic case, the leptokurtic case has only given negative values of B, i.e., 
no other than continuous phase transition occurs. The phase diagram is shown in Fig. [3J on the plane formed by the 
temperature and the PDF width to (in units of J), for some specific values of r > 1. Still, the lines represent numerical 
solutions of Eq. (|14| . while at the same time the points were analytically obtained through a zero-temperature analysis, 
which is going to be discussed shortly. As we have perceived in t he p latykurtic case, the ferromagnetic phase is reduced 
by augmenting r. Similar behavior was found in the Gaussian [17| and the double-Gaussian RFIM [2l| by increasing 
the standard deviation of such PDFs. However, a chief difference emerges. For distributions with fatter tails than the 
Cauchy-Lorentz PDF, the concavity of the critical frontier changes in the high-temperature region. So far as we are 
aware, this is the first time that such a change is observed in the mean-field RFIM phase diagram. 



B. Leptokurtic case: r > 1 
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IV. ZERO TEMPERATURE ANALYSIS 



Moving forwards, we now consider the phase diagram of the model at zero temperature. As in the finite-temperature 
case, we evolve twofold: the platykurtic case and the leptokurtic case, r < 1 and r > 1, respectively. 



A. Platykurtic case: r < 1 

In the limit T — > 0, the free energy and magnetization become 1 , respectively, 
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where 2Fi[-, ■] is the Gauss hypergeometric function [341 ]. In the same way as in the finite-temperature analysis, 
we expand the above magnetization (|20p in powers of m, , so that 
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The continuous critical frontier at zero temperature is obtained for a\ = 1, 
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providing that 6^ < 0, which occurs for all r < 1. The last-mentioned equation allows determining the exact point 
at which the critical frontiers obtained in section 3. A reach the zero-temperature axis (the circles in Fig. [2]). The 
zero-temperature phase diagram is shown in Fig. |4] 



1 For the purpose of obtaining the following expressions we made use of the integrals presented in Ref. [33 



B. Leptokurtic case: r > 1 



In this regime, PDF presents a distinct behavior for 1 < r < 5/3 and t > 5/3. Explicitly, the former case 
corresponds to the case in which the standard deviation is finite and the latter to the case for which the distribution 
has the same asymptotic behavior as the Levy distribution. 



1. Finite standard deviation: 1 < r < 5/3 
For this range of r, the free energy and the magnetization become, 
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respectively. Similarly to the r < 1 analysis, we can expand the magnetization to s , Eq. (|27|l . in powers of : 



m s — a s m s + b s m 3 s + c s ml + 0(m s ), 
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The continuous critical frontier at zero temperature is obtained for a s = 1, 
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as long as b s < 0. In the range 1 < r < 5/3, we notice that the coefficient b s is always negative, indicating the 
occurrence of continuous phase transitions for all values of r. This expression permit us to determine the values of 
a/ J, or equivalently, the values of ui/J [see Eq. (0] at T = of the phase diagrams depicted in Fig. [3l In Fig. IH we 
show the zero-temperature phase diagram, on the plane w/J vs r. 
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FIG. 4: (Color online) Zero-temperature phase diagram separating the Ferromagnetic (F) and the Paramagnetic (P) phases 
for platykurtic (r < 1) and leptokurtic (r > 1) distributions. The horizontal dotted line represents the limiting case r — > -co, 
i.e., the uniform distribution (ui/J = 1), the dashed vertical line represents the limit for finite standard deviation (r = 5/3) 
and the dot-dashed line the limit for finite average (r = 2). The points emphasize the intersection between vertical lines and 
the critical line. They correspond to values of ui/J equal to ^j2/n, \/\/2 and 2 /n, respectively. 

2. Finite width: 5/3 < r < 3 

Mark that in this range we must use to. Thus, as previously, the free energy and magnetization are respectively, 
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Analogously to the above cases, we expand the magnetization m s , Eq. (134p . in powers of m s , 

to s = a s m s + b s m 3 s + c s m 5 s + 0(m r s ), 

with the coefficients, 
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FIG. 5: (Color online) Tri-dimensional phase diagram of the model in the axis temperature, r and ui/J, separating the 
Ferromagnetic (F) and the Paramagnetic (P) phases. We have used a darker color to represent the regime of r (r > 2) in 
which we have determined a divergent free energy per spin analytically found at T — 0. 
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Thus, the continuous critical frontier is given by 
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and we have again verified that b s < for all values 5/3 < r < 3. We can see in Fig. |4]the zero-temperature phase 
diagram in the plane containing the width to and the generalized degree of freedom r. 

In this case, it is worth noticing an important result. In the free energy per spin (|33|) . the integrals are finite only 
for r < 2, i.e., the free energy at temperature equal to zero is not finite for probability density functions broader than 
the Cauchy-Lorentz. Although we do not have an unequivocal physical account for this phenomenon, we introduce 
some insight into this result with the help of the statistical meaning of our distributions. As mentioned in section 2, 
for r > 1, distribution J3J| is understood as a generalization of Student-i for real degrees of freedom according to the 
relation, 

(40) 



1 + m 

The Cauchy-Lorentz distribution, Eq. |3]) with r = 2, corresponds to the case for which the distribution presents 
a divergence in the average but a null average value of the corresponding escort distribution. The divergence of the 
mean value of the free energy for r > 2 emerges from that feature of the property of Eq. Moreover, this divergence 
was experimentally observed in organic charge-transfer compounds [24j |. 

In order to summarize the results presented in the manuscript, we show in Fig. [5] a tridimensional phase diagram 
separating the ferromagnetic (F) and the paramagnetic (P) phases defined by the axis temperature (in units of J), 



r and uj (also in unit of J). We observe a contraction of the ferromagnetic phase for increasing values of r. We have 
spotted the above-described change in the concavity of the critical frontier for r > 2, as well as the dwindling of the 
ferromagnetic phase (for increasing values of r) which in limit r — > 3 turns into the point lu = 0. 

V. CONCLUDING REMARKS 

In this work we have investigated the infinite-range-interaction Ising model in the presence of a random magnetic 
field following a family of continuous probability density functions, defined by a parameter r comprising the re- 
distribution, for r < 1, and the Student-i, for r > 1, which have already found their statistical relevance within 
other contexts of disordered systems. Moreover, specific PDFs like the Gaussian (r — > 1), the uniform (r — > — oo) 
and the Cauchy-Lorentz (t = 2) are obtained thereof. Independently of r, we have observed a continuous phase 
transition with the lessening of the ferromagnetic phase in the ^j- vs j plane that corresponds to the region defined 
by < ^j- < 1 and < j < 1 in the uniform case and to the point j = ^j- = for r — > 3. For r > 2, we have 
noted the appearance of an inflexion point for finite j and ^j- that is also associated with a divergence of the free 
energy per spin at null temperature for which we have provided with an explanation based on the statistical nature 
of distributions that are fatter than the Cauchy-Lorentz. 

As an extension of this work, a numerical study by means of Monte Carlo simulations of the model defined by 
Eqs. ([T]), and ^ in the case of nearest-neighbors interactions is though to bring a better understanding of the 
physical properties of the Ising model in the presence of random magnetic fields that follow continuous probability 
distributions [351 ] . 
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